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We study an extended Shastry-Sutherland model for 
SrCu2 (603)2 and analyze the low lying parts of the energy 
spectrum by means of a perturbative unitary transformation 
based on flow equations. The derivation of the 1-magnon dis- 
persion (elementary triplets) is discussed. Additionally, we 
give a quantitative description (symmetries and energies) of 
bound states made from two elementary triplets. Our high or- 
der results allow to fix the model parameters for SrCu2 (603)2 
precisely: Ji = 6.16(10)meV, x ■- J2/J1 = 0.603(3), Jx = 
1.3(2)meV. To our knowledge this is the first quantitative 
treatment of bound states in a true 2d model. 

PACS numbers: 75.40. Gb, 75.30.Kz, 75.50.Ee, 75.10.Jm 



Presently low dimensional quantum antiferromagnets 
are investigated intensively both in experiment and the- 
ory. Systems that do not show a long ranged ordered 
ground state, so called spin liquids, are particularly in- 
teresting. Besides low spin and low coordination number 
spin liquids are more likely to form in strongly frustrated 
geometries. The recently synthesized antiferromagnetic 
material SrCu2 (303)2 Wpl is a nice two dimensional ex- 
ample since it is a realization of the Shastry-Sutherland 
model |§ (cf. Fig. |i|). 




FIG. 1. A part of the extended Shastry-Sutherland model. The 
circles are S = 1/2 entities. For J3 = this model reduces to the 
original Shastry-Sutherland model. The coupling J\ is assumed to 
be antiferromagnetic. The starting point of our analysis is the limit 
of strong dimcrization (J2, J3 — ► 0). 

In this model frustration plays an essential role. Each 
spin is coupled to pairs of spins on dimcrs. If these pairs 
form singlets the couplings between dimers is without ef- 
fect and the singlet-on-dimers product state is always an 
eigen state and for certain parameters the ground state 



glljjj^. In this dimer phase the system is gaped. The 
Hamiltonian reads 



H^JlJ2S^■ S, +J2 J2 S^ ■ Sk +J3YI S^ ' Si, (1) 
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where the bonds corresponding to interactions Ji , J2 and 
J3 are shown in Fig. H^. 

For larger ratios x = J2/J1 and y = J3/J1 the sys- 
tem transits into a Neel phase, competing with a helical 
phase for certain values of y. Finally there exists a fer- 
romagnetic phase in the region of larger negative x and 
y. A profound analysis of these scenarios can be found in 
Ref. [ 7| . Along the x-axis Koga and Kawakami [|| found 
an additional plaquette-singlet phase, which we see as 
beeing related to the helical phase. 

In this article we start from the dimer phase where an 
elementary excitation is given by breaking up one singlet 
and substituting a triplet instead, which acquires a dis- 
persion by hopping from dimer to dimer. For the next 
higher excitations we focus on bound states formed from 
pairs of the elementary triplets. We will perturbatively 
calculate the low lying excitations of the model about the 
limit of strong dimerization. The problem can be stated 



H 

j7 



= Ho + xHs , with Hs = if 1 



-H, 



(2) 



In the limit of isolated dimers [x ~ 0, with yjx finite) 
H is bounded from below and has an equidistant energy 
spectrum, since Hq simply counts the number of excited 
dimers (up to a trivial constant). Furthermore iJs can 
be decomposed 



Hs = r_i + To + Ti 



(3) 



where Ti creates i (destroys for i < 0) elementary triplets 
(energy quanta). These properties of Hq and H^ al- 
low us to use the perturbative unitary transformation 
H based on flow equations [nO|. This technique enables 
us to link smoothly and uniquely if at x ^ to an ef- 
fective iieff conserving the number of triplets on dimers 
[ifcffjifo] — 0. This permits a clear distinction between 
the ground state (without triplets), the 1-triplet sector, 
the 2-triplet sector etc.. Thus the effective Hamiltonian 
is block diagonal. The explicit form is given by 



Hoff = 7^0 + 51^'' H 



C{m)T{m) 



(4) 



fe=l \m\=k,M{m)=0 



where m is a vector of dimension k of which the compo- 
nents are in {±1, 0}; M{m) = signifies that the sum of 
the components vanishes which reflects the conservation 
of the number of energy quanta. The operator product 
T{rn) is defined by T™ = Tm-i^Tm2 • • 'Tm^- The coeffi- 
cients C{rn) are generaUy vafid fractions, which we com- 
puted up to order k — 15. 

In terms of HeS it is easy to show that a hopping, 
displacing an elementary triplet by a finite distance, 
starts in sixth order in x, while in the 2-triplet sector 
correlated hopping starts in second order (see uB and 
Refs. therein), explaining the rather fiat 1-triplet dis- 
persion in contrast to the much stronger pronounced 2- 
triplet dispersion as found by Kageyama et al. by INS 
measurements 1 12 . By means of degenerate perturbation 



theory Momoi and Totsuka |]13[ also derive an effective 
Hamiltonian for the Shastry-Sutherland model. Their 
third order result also shows the significance of correlated 
hopping in the multi-triplet dynamics of the model. 
Let |r) = |ri, r2) denote the state of the system with one 
triplet on the dimer at r and singlets on all other sites. 
The amplitude t°;_^ for a triplet-hopping from site r to 
site r' is then given by 



o(r) 



= (r'|Heff|r) 



(5) 



where the upper index o(r) G {v, h} allows to distinguish 
whether the hopping started on a vertically oriented (v) 
or a horizontally oriented dimer (h). Further we choose 
to split the hopping amplitudes into a net part tg and a 
deviation part dig (s = r' — r) 



o(r) _ 



U + e'^'dU 



with Q = (tt, tt). 

Since Hcs conserves the number of triplets one has 



ifeff|r) = $:C^--)|r + r') 



We introduce the Fourier transformed states 



\a,k) 



m^^ 



J{k+aQ)r 



(6) 



(7) 



(8) 



with the total number of dimers N, the new quantum 
number a € {0, 1} reflecting the sub lattice structure and 
k a vector in the magnetic Brillouin zone (MBZ). Hcs 
acts as a 2 X 2 matrix on the states \a, k) and |1 — cr, k). 
Its diagonalization yields 



wi/2(k) 



flo + ai 1 



2 

a)o(k) 

Here we have defined 



±-V(ao-ai)2 + 462. (9) 



to + 2^tr cos((k + uQ)r) 

r>0 

& = 2 2_\ dtr cos(kr) , 



and 



(10) 



r>0 

r\-\-r2 even 



with r > if ri > 0, or ri = but r2 > 0. Thus the 
1-triplet dispersion splits into two branches. We want to 
point out, however, that at k = and on the borders of 
the MBZ the two branches fall onto each other leading to 
a 2-fold degenerate dispersion at these points. This can 
be derived by showing that the square root in Eq. (|9|) 
vanishes at these points ||lj]. An analogous degeneracy 
is noticed in the 2-triplet sector. In Ref. ||ll| , [l4| we give 
a detailed analysis of the symmetries of the model (2d 
space group p4mm with underlying point group 4mm) 
and show that the degeneracies are due to glide line sym- 
metries. 

We calculated the amplitudes t^ and dtr (and therefore 
the dispersion) as exact polynomials in Ji , a: and y up to 
and including 15**^ order. 

Expanding the square root in Eq. (||) about the limit of 
vanishing x and y produces terms oc x'^y^ with a -\- (3 > 
10. Hence the energy splitting starts in 10**^ order only. 
It is negligible for all reasonable values of x and y. 
By substituting y = in wo(k) we verify the decimal 
numbers Weihong et al. pq ] obtained previously. 
For three different sets of Ji , x and y the 1-triplet disper- 
sion is plotted in the inset of Fig. |. ESR @, FIR |l|] 
and INS Q data suggest a value of u;(0) = 2.98meV. At 
finite k we have to rely on the INS measurement, which 
have rather large errors. We get a very good agreement. 
But it is not possible to fix the model parameters unam- 
biguously from information in the 1-triplet sector alone. 
Thus we continue our analysis in the 2-triplet sector. 
The dynamics of two triplets at large distances is gov- 
erned by 1-triplet hopping. At smaller distances a 2- 
particle interaction occurs additionally given by WOi:d;r,d' 
(W^t>:d;r.d') Starting with one triplet on a horizontal (ver- 
tical) dimer and another at distance d. The action of 
ifeff is to shift the triplets to r and to r -|- d'. Nothing 
else is possible due to triplet number conservation. Since 
the total spin is conserved (5* S {0, 1,2}) the distances 
are restricted to d,d' > 0, because the exchange parity 
is fixed. 

The coefficients W for S* € {0, 1, 2} and the 1-triplet hop- 
ping yield the complete 2-particle dynamics. We compute 
W up to x^"^ , the coefficients for the lowest- lying states 
even up to x^"*. 

Analogously to the 1-triplet problem we use the following 
basis for the 2-triplet states 

|k, d, a) := iV-i/2 ^ e(*(k+.Q)(r+d/2)) j^.^ ^ ^ ^^ ^ (^^^ 

r 

where k is the conserved total momentum in the MBZ 
applying due to the two sub lattices; |r,r + d) denotes 



the state with triplets at r and at r + d (d > 0). 
Since 1-triplet hopping is of higher order than the in- 
teraction an analytic expansion for the energies of the 
bound states is possible. At finite order in x only config- 
urations contribute where the two triplets are not too 
far away from each other. Of course, higher orders 
imply larger, but still finite distances. In particular, 
the energies of the four states which evolve from neigh- 
boring triplets can be computed very well since their 
leading interaction is linear. Investigating the matrix 
elements shows that it is sufficient to study the dis- 
tances d e {(0, 1),(1,0),(1,±1)} for order 5. To x'^'^ 
only d e {(1, ±2), (2, ±1), (0,2), (2,0), (2, ±2)} must be 
added. So, for given total momentum only a finite 8x8 
or 24 X 24 matrix has to be analyzed. Furthermore, the 
elements connecting shorter distances to longer distances 
and the elements among longer distances do not need to 
be known to very high order. 

We have analyzed the dispersions in x^ of the four 
states bound linearly in x in the MBZ. Fukumoto's 
results are mostly confirmed |18|. The dispersion of 
bound states starts only in x^ (contrary to x^ claimed 
in Ref. (l^). At particular points of high point group 
symmetry ((0,0),(0,7r)) the Hamiltonian splits into six 
blocks corresponding to different representations of the 
square point group 4mm. At these points the analy- 
sis up to x^'^ is carried out. The symmetries are clas- 
sified according to the irreducible representations (four 
ID, one 2D) of the point group 4mm ri(l),r2(a;^ — 
y'^),'^3{xy),r4{xy{x'^ -y'^)),T5{x,y) where simple poly- 
nomials are given in brackets to show the transformation 
behavior. 
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FIG. 2. Energy of the lowest-lying 5 = states. Curves refer to 
k = except the dashed-dotted one. The dotted curve displays the 
continuum at 2 A. Inset; 1-triplet dispersion for various values of 
(Ji/meV;x;y): (6.56;0.615;0) dashed-dotted line, (6.67;0.59;-0.05) 
dashed line, (6.16;0.603;0) solid line. 

The extrapolated energies are depicted in Figs. S (5* = 
0) and H (5 = 1) as functions of x. The double de- 



generacy for k = (0, 7r) results from the same symmetry 
reasons as described for the 1-magnon dispersion. The 
dashed-dotted curve at (0, tt) has to be compared to the 
solid and the long-dashed curve to assess the dispersion 
of these two modes from to (0,7r). While for 5' = 
this dispersion always has the expected behavior with 
u!{0) < ix'((0,7r)) the energies for S = 1 are reversed for 
small values of x (cf. p8|). Only above x ~ 0.55 the 
relation a;(0) < w((0,7r)) holds for S = 1. 
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FIG. 3. As in Fig. Q for 5 = 1. Inset: Magnetic susceptibility. 
Theory (dashed) with directional rms average g = 2.13 |16[ , x,Ji 
as in Fig. eI experiment (solid) on powder M. 

For 5" = 0, the lowest mode vanishes at the same 
X as does the elementary triplet gap A. So no addi- 
tional instability occurs for S — 0. This provides evi- 
dence against the competing singlet phase as presumed 
in Ref. [gj. There is, however, a salient instability for 
S = 1 (Fig. ^) at x = 0.63. This comes as a surprise 
since one expects in antiferromagnets binding effects to 
be largest for S = (discussed in more detail in Ref. |Q ) . 
Assuming A — 2.98meV as above and uj\s=i = 4.7meV 
from ESR |§, FIR and INS |l|l we can determine 
the model parameters for SrCu2(B03)2 precisely {y = 0): 
X = 0.603(3) and Ji = 6.16(10)meV. 
The inset in Fig. g shows that the 1-magnon dispersion 
agrees very nicely with experiment for these values. Fur- 
ther, the energy of the Fa singlet matches the 30cm~"'^ 
peak in Raman scattering |19] perfectly. The r4 singlet at 
25cm~^ is forbidden by symmetry, since the Raman op- 
erator is Fa at T=0 ||ll|. Calculating the next F3, 5* = 
bound state (not shown) yields 45cm~^ in good agree- 
ment with the experimental 46cm~^ line, too. 
We conclude that the 2D model (Fig. 1^) explains the 
low-lying excitations of SrCu2 (603)2 perfectly. Ther- 
modynamic quantities like the susceptibility x(T') require 
the inclusion of the interplane coupling J± which is fully 
frustrated not changing the dimer spins |Q. The 3D 
X3D is computed from X2B on the mean-field level as 



XsD = X2D + 4>/_L (inset in Fig. ||). Fits to the experi- 
ment Ig] give J_L = 1.3(2)meV leading to a curve agreeing 
without flaw above 40K. Our value for J± is significantly 
higher than the one in Ref. [^ due to different values of 
X and Ji. 

Summarizing, we presented the first quantitative descrip- 
tion of 2-particle bound states in 2D. 
An unexpected instability for the S = 1 2-triplct bound 
state is found at x w 0.63 indicating a transition to a 
triplet condensate probably related to the helical phase 
found previously |21 ^. We conjecture that this transi- 



tion is first order occuring at lower x than assumed so 

far. 

The symmetries of experimentally relevant states were 

determined. The reliability of the high order results 

allows to fix the experimental coupling constants very 

precisely (Ji = 6.16(10)mev, J2/J1 = 0.603(3), J_l = 

1.3(2)meV). Thereby, various experiments (ESR, FIR, 

INS, Raman, x(T)) are explained consistently. 

The work is supported by the DFG in SFB 341 and in 

SP 1073. 
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